home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dtbcon.z / dtbcon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDTTTTBBBBCCCCOOOONNNN((((3333FFFF))))                                                          DDDDTTTTBBBBCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DTBCON - estimate the reciprocal of the condition number of a triangular
  10.      band matrix A, in either the 1-norm or the infinity-norm
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DTBCON( NORM, UPLO, DIAG, N, KD, AB, LDAB, RCOND, WORK, IWORK,
  14.                         INFO )
  15.  
  16.          CHARACTER      DIAG, NORM, UPLO
  17.  
  18.          INTEGER        INFO, KD, LDAB, N
  19.  
  20.          DOUBLE         PRECISION RCOND
  21.  
  22.          INTEGER        IWORK( * )
  23.  
  24.          DOUBLE         PRECISION AB( LDAB, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      DTBCON estimates the reciprocal of the condition number of a triangular
  28.      band matrix A, in either the 1-norm or the infinity-norm.
  29.  
  30.      The norm of A is computed and an estimate is obtained for norm(inv(A)),
  31.      then the reciprocal of the condition number is computed as
  32.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  33.  
  34.  
  35. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  36.      NORM    (input) CHARACTER*1
  37.              Specifies whether the 1-norm condition number or the infinity-
  38.              norm condition number is required:
  39.              = '1' or 'O':  1-norm;
  40.              = 'I':         Infinity-norm.
  41.  
  42.      UPLO    (input) CHARACTER*1
  43.              = 'U':  A is upper triangular;
  44.              = 'L':  A is lower triangular.
  45.  
  46.      DIAG    (input) CHARACTER*1
  47.              = 'N':  A is non-unit triangular;
  48.              = 'U':  A is unit triangular.
  49.  
  50.      N       (input) INTEGER
  51.              The order of the matrix A.  N >= 0.
  52.  
  53.      KD      (input) INTEGER
  54.              The number of superdiagonals or subdiagonals of the triangular
  55.              band matrix A.  KD >= 0.
  56.  
  57.      AB      (input) DOUBLE PRECISION array, dimension (LDAB,N)
  58.              The upper or lower triangular band matrix A, stored in the first
  59.              kd+1 rows of the array. The j-th column of A is stored in the j-
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDTTTTBBBBCCCCOOOONNNN((((3333FFFF))))                                                          DDDDTTTTBBBBCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              th column of the array AB as follows:  if UPLO = 'U', AB(kd+1+i-
  75.              j,j) = A(i,j) for max(1,j-kd)<=i<=j; if UPLO = 'L', AB(1+i-j,j)
  76.              = A(i,j) for j<=i<=min(n,j+kd).  If DIAG = 'U', the diagonal
  77.              elements of A are not referenced and are assumed to be 1.
  78.  
  79.      LDAB    (input) INTEGER
  80.              The leading dimension of the array AB.  LDAB >= KD+1.
  81.  
  82.      RCOND   (output) DOUBLE PRECISION
  83.              The reciprocal of the condition number of the matrix A, computed
  84.              as RCOND = 1/(norm(A) * norm(inv(A))).
  85.  
  86.      WORK    (workspace) DOUBLE PRECISION array, dimension (3*N)
  87.  
  88.      IWORK   (workspace) INTEGER array, dimension (N)
  89.  
  90.      INFO    (output) INTEGER
  91.              = 0:  successful exit
  92.              < 0:  if INFO = -i, the i-th argument had an illegal value
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.